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We consider the number of configurations of a surface in two dimensions that has a prescribed 
length and encloses a prescribed perimeter with respect to a baseline. An approximate analytical 
treatment in a semi-continuum compares favourably with results from an exact algorithm for the 
discrete lattice. This work is relevant for finding the entropy associated with macroscopic configu- 
rations of such systems as domain growth problems, evaporation-deposition problems, membrane 
physics, or polymer physics. 
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I. INTRODUCTION 



II. DEFINITION OF THE PROBLEM 



In statistical mechanics, one very often performs aver- 
ages over ensembles of microscopic configurations. Since 
this is usually very difficult, it is common to replace the 
ensemble of microscopic quantities by an ensemble of 
macroscopic quantities. For instance, a spin Hamiltonian 
may be replaced by one that depends only on magnetiza- 
tion. This paper, we addresses the problem of a surface 
in two dimensions with fixed endpoints where the macro- 
scopic quantities are the perimeter of the surface and the 
area enclosed by it. Such a problem would arise, for ex- 
ample, if one is given a Hamiltonian for a system that 
consists of a area term and a surface term, i.e. 



n= Jfs{S)-PfaiA) 



(1) 



where S is the perimeter of the surface and A is the area 
between it and some substrate, J is the surface tension, 
P the pressure, and fs, fa general functions of S and 
A respectively. These types of Hamiltonians could arise 
in domain roughening problems wetting problems, 

evaporation-deposition problems, the physics of mem- 
branes, or polymer physics [||J^. We were motivated by 
an investigation of the random field Ising model |7 . 

To find the free energy of such a system and thus 
any relevant macroscopic quantities at a finite temper- 
ature, one would like to know the associated entropy 
S = fcs lnfiL(S', A). Here, 1)^(5, A) is the number of 
surface configurations with the given perimeter S and 
enclosed area S and which begin and end at the ends of 
a line of length L; this is the quantity that shall be of 
interest to us. 

A slightly different version of this problein|^ has been 
solved exactly by Owczark and Prellberg |^ by finding 
the generating function of fi. However, the solution was 
in terms of the q-series deformation of a hypergeometric 
function. While informative, this result is very complex 
and difficult to use in our analysis. Our aim is to find a 
simple form for that can be easily manipulated. 

We use an approximation to translate the problem into 
one of simple Riemannian geometry that can be solved 
using the familiar tools of vector analysis. The result is 
the pleasingly elementary form 



n{S,A)^f[g'S^ 



A^ 

£2 



(2) 



where /, C and g are functions of L only. This result 
compares favourably with exact numerical results. 



^They studied "bar graph polygons" wherein the domain 
wall heights are not allowed to go below the baseline. In our 
model, the heights of the surface can be both positive and 
negative. 



We consider a surface with fixed endpoints such that 
it has an end to end distance of L. The surface may take 
any path between these endpoints as long as there are no 
overlaps or overhangs {i.e. it is a directed random walk). 
This should be a reasonable approximation in the limit 
of weak disorder. We label the straight line between the 
endpoints with the integer i S [0, L] and the height of the 
surface at each point hy hi, a single valued real number. 
Figure |^ illustrates these definitions. The quantities A 
and S may be written as functions of these variables. 



A = Y,h^ 

i=0 
L 

S ^^^\hi - hi-i\ 



(3a) 
(3b) 



where S is the actual length of the surface less L (since 
any surface that covers the baseline distance must be at 
least L long). 

The problem is simplified somewhat by the introduc- 
tion of the variables Ri — hi — Like hi, Ri is a real 
number with i S [1,L] an integer. Under this change of 
variables, equation (||) becomes 



A = J2iL~i)R^ 

1=1 

L 

s = j2m. 



(4a) 



(4b) 



i=l 



There is also an additional constraint that the surface 
finish at height h^ = 0. This introduces the auxiliary 
equation 



^ i?, = 0. 



(5) 



Consider an ensemble of surfaces all with various {A, S) 
and subject only to the constraint (||). In terms of the 
definitions above, we wish to know how many different 
values of Ri satisfy equations ( ^ ) , ( |4b| ) , and (^) for fixed 
A and S. We solve this problem in the next section. 



III. SOLUTION 

The solution is further complicated by the presence of 
the absolute value function in (4b) and so we will make 
a further approximation to simplify matters. It is hoped 
that this will not change the essential character of the 
solution. Thus, we replace equation (^) by 



g{LrS' = J2LR^ 



(6) 
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This completes the approximations in the solution. That 
the form factor g should only be a function of L follows 
from a simple scalingargument. 

Clearly, equation (|6|) implies that the vectors Ri lie on 
a. L — 1 dimensional hypersphere embedded in L dimen- 
sions. The other two equations and (||) introduce 
hyperplanes that pass through the hypersphere. Thus, 
vectors that satisfy all three equations lie on the inter- 
section of the hypersphere with these two hyperplanes. 
This intersection is itself a hypersphere, but of dimen- 
sion L — 3. Thus, the problem is reduced to finding the 
measure of the vectors which lie on the surface of this 
hypersphere. The details of this calculation are included 
in appendix ^ and the result is 



nL(5,A)=/(max(g^52_^^0 



(7) 



L-2 



!(L-2) 



It is worth noting that the problem posed by equa- 
tions (H), (^a]), and ^ with g{L) = 1 is in itself worthy 
of consideration, and may relate to other models in sta- 
tistical physics. However, since we were motivated by the 
random field Ising model wherein ( ^ ) arises, we set g{L) 
so that maps as closely as possible back to (|^. We 
describe how to do this in the following section. 



IV. THE FORM FACTOR g{L) 

One can employ simple scaling arguments on to 
deduce that g{L) ~ cL~^l'^ for large values of L, where 
we would expect c to be of order unity. But, a better 
estimate of g can be made. Consider the set of vectors 
that satisfy (|^), that is, the vectors on an L dimensional 
hypersphere of radius r ~ g{L)A. To satisiy (|^) we 
must set A such that — ^ ^'^^ vectors on 

this hypersphere. Obviously, we cannot satisfy this crite- 
rion for all vectors, so we demand that it be satisfied on 
average. This presumes, of course, that all the vectors 
on the hypersphere have equal importance in O^, which 
is a reasonable first approximation. Thus, we set 



(El 



= A 



(8) 



where (Q) denotes the average value of a quantity Q over 
the hyperball of radius r. It is a straightforward piece of 
mathematics to perform this calculation, the details of 
which are included in appendix H. The result is 



'L±l^ 

- 2 } 



(L/2)! (1 -f 1/L) , 



(9) 



Using Stirling's approximation, we find that for large L 
this tends asymptotically to 



n 

2L' 



(10) 



which agrees with our earlier expectations. 



V. THE NUMERICAL TREATMENT 

In order to test (and possibly improve upon) the an- 
alytic results obtained above, an algorithm was devised 
to scan over all the random walks with perimeter S less 
than some preset maximum A*. To understand this al- 
gorithm, it is helpful to think of S as the length of a 
piece of string that we have to play with. This piece of 
string must start and finish on the baseline but may do 
anything in between (subject to the restriction that there 
are no overhangs or overlaps). But, the amount of string 
we have used up in our previous steps determines how 
much we have left to play with in the following steps. 
For instance, at the first step, the string may have any 
height Ro G [—A*/2,A*/2]. In general, if we are at posi- 
tion i and at height hi = X]]=i j ^^'^ have already used 
up a length of string Ai — l-^j I; then this limits the 

range of values Ri may have. Namely, 



R^ e 



'A* + Ai A* - Ai 



(11) 



Figure |2| illustrates this. 

Furthermore, we only need to scan over various values 
of Ri for i < L~l since R^ is determined by the condition 
— 0. A general outline of the algorithm is 
given below: 

1. Set Ri to -A*/2 

2. Set all Ri for i = 2, . . . , L — 1 to their minimum 
possible value and determine Rl 

3. Measure S and A. Increase il{S,A) by one. 

4. Increase Rl-i by one. If exceeds its maximum 
value then increase Rl-2 by one and set Rl-i to 
its minimum value. If this causes Rl-2 to exceed its 
maximum value, then increase Rl-3 by one and set 
Rl-2 to its minimum value. Et cetera. 

5. Repeat steps 3 and 4 until i?i = A* /2 

A recursive program was written in C to implement 
this algorithm, and the program was run on Sun work- 
stations. Data was obtained for A* as large as 60 and 
for L as large as 10. It becomes increasingly difficult to 
get data for larger L since the number of random walks 
seems to grow geometrically with L. Indeed, we predict 
that 



n{S,A)dV {L^2)^A^a^ 



where a is a constant of order unity. 



(12) 
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VI. COMPARISON OF NUMERICAL AND 
ANALYTIC RESULTS 

When viewed side by side, the two results seem to look 
very similar (see Figure |3|). However, the scales do not 
agree. This is attributed to the fact that the numerical 
routine investigates only discrete values of Ri whereas in 
the analytic treatment Ri was allowed to be continuous. 
In both treatments L is divided into discrete intervals, so 
this is not a cause for worry (as long as we have L > 3). 
From the point of view of our hypersphere argument, the 
analytic treatment has given the measure of all vectors on 
the surface, whereas if the components of the vectors are 
constrained to be integers, there are fewer vectors close 
to the surface of the hypersphere. Nowhere is the effect 
of the discretization more noticeable than in the fact that 
il{S,A) = on a discrete lattice for odd S. This simply 
says that we need as many steps up as down if we are to 
return to the origin. We do not have this same constraint 
in the semi-continuous model. Nevertheless, the results 
still compare very favourably on various levels (up to this 
scaling factor). Three quantities were chosen for the pur- 
poses of comparison, which we list below together with 
their predicted values: 



The maximum value of Q for a given S and L, 



Mm^nL{A,0)^ f{L) {g{L)A) 



L-3 



(13) 



• The zeroth moment of the distribution with respect 
to A, 

Mo = J n{S, A)dV = Cf{L) {9{L)Af-^ T{L) 

(14) 

where 

22i-5(x,_3)!(X-3)! 



T(L) 



(2L-5)(2L-6)! 



(15) 



• The second moment of the distribution with respect 
to A, 



M2 = y niS,A)A^dV ^C^'fiL) {g{L)AY 



T(£) 

2L-3' 
(16) 



Given that we expect these forms to arise, we generate 
these quantities from our data and fit them to a power 
law behaviour in S for various values of L. The exponent 
and coefficient of this fit are then compared to the ana- 
lytic predictions listed above. The results are shown in 
Figure where one can see that the exponents, but not 
the coefficients, agree well (as expected from the above 
discussion). However, the coefficients do seem to have 
the correct scaling behaviour in L. 



VII. CONCLUSIONS 

We have investigated the question of the entropy of 
a two dimensional surface with fixed volume and area 
from both the numerical and analytic viewpoints. From 
the analytic side, we found an elegant expression for the 
semi-continuum by using a single approximation. From 
the numerical side, we generated exact results that com- 
pared favourably with the analytic expression. However, 
there were certain differences noted between the two and 
these were attributed to differences between the discrete 
and semi-continuous formulations of the problem. 

Nevertheless, the solution captures the important char- 
acteristics of the numerical work. It is hoped that this 
work will be of value to finite temperature formulations 
of problems in many areas of statistical mechanics where 
the question may be posed in terms of such macroscopic 
quantities as the length of a curve and the area enclosed 
by it. 
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APPENDIX A: CALCULATION OF Q.{A, V) 

We wish to find the measure of all vectors that sat- 
isfy equations (^), (||), and (|^). To simplify matters, let 
us introduce the dimensionless variables r[ — Ri/{gA). 
Working with these vectors instead of Ri introduces a 
change of measure and since the object we are looking at 
is clearly a L — 3 dimensional hypersphere, this change 
of measure is J = {gAV^^. Written in these new coor- 
dinates the equations (pa|), (H), and (^) become 



'2 1 

r, = 1, 



E 

i=l 

L 

J2{L- 1 + 1)4 ^v, 



(Ala) 
(Alb) 
(Ale) 



where v = V/g{L)A. 

This change of measure associated with the change of 
variables is associated with the scaling behaviour of the 
equations. We can make this more explicit by writing 



nL{A,V) = mL{l/g{L),V/g{L)A). 



(A2) 
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Similarly, we could have written 

nL{A, v) = v'^-^nL{A/v,i). 



(A3) 



It is important that our final solution obey these two 
relations. 

The path to th e so lution is made clearer if we rewrite 
the equations in ( |Al[ ) so that they more closely resemble 
the vector fomalism. To achieve this, let us introduce the 
vectors 



7 '2i ■ 



^1} 



1 ={1,1,...,!} 

i' = {i-l,L-2, 



,0} 



so that (Al) becomes 



¥\\ = i 

r • i' = 

r' -L' ^ WL'W cose 



(A4a) 
(A4b) 
(A4c) 

(A5a) 
(A5b) 
(A5c) 



where \\LT ^ L(i - l)(2i - l)/6. 

Figure |^ gives a pic torial representation of these equa- 



tions. Equation ( A5a) specifies that the vectors lie on a 



unit hypersphere, 5", of dimension L—1. Equation (A5b) 
is an orthogonality condition and is used to lower the di- 
mension of the problem by one. Finally, equation (A5c) 



specifies that the vectors are at an angle 9 to the vector 
L'. 

We will now employ a coordinate rotation and project 
our problem onto the L—1 dimensional space othogonal 
to 1. The rotation is chosen to simplify the projection, 
e.g. we let the z axis lie parallel to the vector 1 so that 
r' (r,0) and L' 

Now f and L are both L—1 dimensional vectors and 
Lz is the component of L' parallel to 1. 



L' -i _L{L- 1) 
||f|| " 2^/L 



(A6) 



Equation (A5c) retains its familiar form in this new, un- 
primed, coordinate system 



' ■ L = \\L\\ COS0 = V. 



(A7) 



We want a measure of the set of all unit vectors f that 
satisfy this dot product. These vectors form a hypercone 
with angle of opening 29 and sides of unit length. An 
appropriate measure is the area of the L — 2 dimensional 
hypersphere that the tips of these vectors trace out. In 
general, the area of a d dimensional hypersphere of radius 
r is 



,d-l 



(d/2)! • 

Letting d L — 2 and r — > sin 9 we arrive at 



(A8) 



n{i,v) = TT- 



L-2 



I (i - 2) (1-cos^ 61)^. (A9) 



We substitute for the value of cos 9 using equation ( [A7| ) 
and simplify the notation by letting C = \\L\\, i.e. 



C'^\\LT-Ll 

_ L{L-l)(2L-l) L{L-lf 
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cos 9 = —. 
C 



(AlO) 
(All) 



Using the scaling relation ( A2) we arrive at our final form 



L-2 



£2 



UL-2) 



(A12) 



It is simple to check that this form satisfies the scaling 
relations ( [A2| ) and (A3). Obviously, for physical reasons, 
we set n{A, V) = OilV> g[L)AC. 



APPENDIX B: CALCULATION OF g{L) 



Following the philosophy outlined in section IV we wish 



to find the sum of the absolute values of the components 
of a vector that lies on an L dimensional hyphersphere of 
radius r. Written in terms of spherical polar coordinates, 
these components are 

xi = r cos 6*1 

X2 = r sin 9i cos 02 



i-l 

Xi = rY\_ sin 9j cos ( 



XL 



L-1 

r Y\ sin 9j 



(Bl) 



We wish to consider (X^iLi \xi\) but we can take ad- 
vantage of the spherical symmetry to write this as L{xi). 
All that remains is to tackle the integral 



G<19iA92 ...d9L- 



L-l 

n 



smf 



(B2) 



where G = r^^^-'^^ Jjf^^^ sm^^^''-^^ 9, is the metric of 
the surface of this hypersphere. Expanding out the inte- 
grals. 
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/ = 2TT / sm9ism^-'-^ OidOi sin9L-iA9L-i. 

(B3) 

The remaining integral may be done by comparison to 
the area integral of the d dimensional unit sphere: 

AM = 2^ n sm'^-'-' e,de, = ^ (B4) 
since 

I = 2l[ r'e,d0, = ^2^^. (B5) 

Thus, the normalized average is (X^^i |a;i|) = 
LI/Ao{L) = A. Putting r = Ag{L) we get the final 
form 



FIG. 4. Plots of the exponent and coefHcient of Mm, Mo 
and M2. The data points are those values found from nu- 
merical work while the solid linos represent the theoretical 
prediction. The error bars represent the statistical error from 
the power law fit. 

FIG. 5. The vector space that 1 and L' live in, showing 
the hyperspheres S and S' . 
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FIG. 1. Directed surface with area S and perimeter S. 



FIG. 2. The drawing on the left shows the maximum value 
Ro can have, thus using up all the string. The drawing on the 
left show the maximum value Ri can have, so that we use up 
all the available string at this step. 



FIG. 3. Sample data set from the numerical routine, com- 
pared with the results of the analytic calculation 
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Figure 1 
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Figure 2 
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